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^[ ', Abstract 
(N 

The aim of these lecture notes is to familiarize graduate students and beginning post- 
graduates with the basic ideas of linear cosmological perturbation theory and of structure 
formation scenarios. We present both the Newtonian and the general relativistic approaches, 
I derive the key equations and then apply them to a number of characteristic cases. The gauge 

Tij- . problem in cosmology and ways to circumvent it are also discussed. We outline the basic 

framework of the baryonic and the non-baryonic structure formation scenarios and point 
. out their strengths and shortcomings. Fundamental concepts, such as the Jeans length. Silk 

damping and coUisionless dissipation, are highlighted and the underlying mathematics are 
presented in a simple and straightforward manner. 

^ ■ 

O ■ 1 Introduction 

^ ' Looking up into the night sky we see structure everywhere. Star clusters, galaxies, galaxy 

clusters, superclusters and voids are evidence that on small and moderate scales, that is up to 
10 Mpc, our universe is very lumpy. As we move to larger and larger scales, however, the universe 
^ ' seems to smooth out. This is evidenced by the isotropy of the x-ray background, the number 

■ counts of radio sources and, of course, by the high isotropy of the Cosmic Microwave Background 

(CMB) radiation. The latter also provides a fossil record of our observable universe when it was 
roughly 10^ years old and about 10^ times smaller than today. So, the universe was very smooth 
at early times and it is very lumpy now. How did this happen? Although the details are still 
elusive, cosmologists believe that the reason is "gravitational instability". Small fluctuations in 
the density of the primeval cosmic fluid that grew gravitationally into the galaxies, the clusters 
and the voids we observe today. The idea of gravitational instability is not new. It was first 
introduced in the early 1900s by Jeans, who showed that a homogeneous and isotropic fluid is 
unstable to small perturbations in its density [Q]. What Jeans demonstrated was that density 
inhomogeneities grow in time when the pressure support is weak compared to the gravitational 
pull. In retrospect, this is not surprising given that gravity is always attractive. As long as 
pressure is negligible, an overdense region will keep accreting material from its surroundings, 
becoming increasingly unstable until it eventually collapses into a gravitationally bound object. 
Jeans, however, applied his analysis to a static Newtonian fluid in an attempt to understand 
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the formation of planets and stars. In modern cosmology we need to account for the expansion 
of the universe as well as for general relativistic effects. 

Despite the lack, as yet, of a detailed scenario, the rather simple idea that the observed 
structure in our universe has resulted from the gravitational amplification of weak primordial 
fluctuations seems to work remarkably well. These small perturbations grew slowly over time 
until they were strong enough to separate from the background expansion, turn around, and 
collapse into gravitationally bound systems like galaxies and galaxy clusters. As long as these 
inhomogeneities are small they can be studied by the linear perturbation theory. A great ad- 
vantage of the linear regime is that the different perturbative modes evolve independently and 
therefore can be treated separately. In this respect, it is natural to divide the analysis of cosmo- 
logical perturbations into two regimes. The early phase, when the perturbation is still outside 
the horizon, and the late time regime when the mode is inside the Hubble radius. In the first case 
microphysical processes, such as pressure effects for example, are negligible and the evolution of 
the perturbation is basically kinematic. After the mode has entered the horizon, however, one 
can no longer disregard microphysics and damping effects. 

In these lectures we will a priori assume the existence of small inhomogeneities at some initial 
time in the early universe. A cosmological model is not complete, however, unless it can also 
produce these seed fluctuations through some viable physical process. Inflation (Guth (1981); 
Linde (1982) |^]) appears to be our best option as it naturally produces a spectrum of scale- 
invariant gaussian perturbations. Topological defects (Kibble (1976) [^), such as cosmic strings, 
global monopoles and textures, offer a radically different paradigm to inflation for structure 
formation purposes. They have fallen out of favor, however, as their observational situation 
looks unpromising. 

Understanding the details of structure formation requires, among other things, knowledge 
of the initial data. Structure formation, or galaxy formation as it is sometimes referred to, 
began effectively with the end of the radiation era at matter-radiation equality (teq — 4.4 x 
10^^^(0/1^)"^ sec, where Q is the density parameter of the universe). Thus, the start of the 
matter era also signals the beginning of structure formation. If we were ever to find out the 
details of how the structure in our universe formed we need to know the initial data at that epoch. 
The necessary information includes: (i) the total amount of the non-relativistic matter; (ii) the 
composition of the universe and the contribution of its various components to the total density, 
namely fib from baryons, from relativistic particles, r^wiMP from relic WIMPs,[| JIa from a 
potential cosmological constant etc; (iii) the spectrum and the type (i.e. adiabatic or isothermal) 
of the primeval density perturbations. Given these one can, in principle, construct a detailed 
scenario of structure formation, which then will be tested against observations. The importance 
of specifying the initial conditions is paramount, since inverting present observations to infer the 
initial data is unfeasible after all the astrophysical filtering that has taken place. Speculation 
on the history of the early universe, backed by recent observations provide some "hints" as to 
the appropriate initial data. They point towards Q = 1 from inflation; r^wiMP — 0.3 (including 
a small baryonic contribution) and ~ 0.7 from nucleosynthesis, inflation and the supernovae 
redshift measurements; and adiabatic fluctuations with a Harrison- Zeldovich spectrum Q from 
inflation. 

The layout of these notes is as follows. In Sec. 2 we present the Newtonian treatment 
of linear cosmological perturbations, discuss issues such as pressure support and the "Jeans 

^Weakly Interacting Massive Particles (WIMPs) are stable non-baryonic species left over from the earliest 
moments of the universe. 
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length" and provide the key results. The general relativistic analysis is outlined in Sec. 3 and 
the basic linear equations are derived. We also give a brief discussion of the "gauge problem" 
and provide characteristic solutions of the relativistic approach. In Sec. 4 we discuss entirely 
baryonic structure formation scenarios, emphasizing the collisional damping of adiabatic density 
perturbations. Non-baryonic "hot" and "cold" dark matter models are presented in Sec. 5, 
together with their advantages and shortcomings. The aim of these lectures notes is to provide 
the basic background to graduate students in physics and astronomy as well as to beginning 
postgraduates. We would like to familiarize the newcomer with fundamental concepts such as 
gravitational instability, the Jeans length, collisional and collisionless damping. The necessary 
mathematics are also provided in simple and straightforward manner. Overall we want to give a 
brief but comprehensive picture of the linear regime, so that the interested student will feel more 
confident when looking at more sophisticated treatments. For further details we refer the reader 
to some of excellent monographs that now exist in the literature (see |8| for a list of them). 
The lectures do not require a particularly specialized background, although some knowledge of 
cosmology and general relativity will be helpful. 

2 Linear Newtonian perturbations 

The Newtonian theory, as a limiting approximation of general relativity, is only applicable to 
scales well within the Hubble radius where the effects of spacetime curvature are negligible. ^ 
Even in this context, however, one can only analyze density perturbations in the non-relativistic 
component. Perturbations in the relativistic matter, at all scales, require the full theory. After 
equipartition the universe is dominated by non-relativistic pressureless matter, which is com- 
monly referred to as "dust". Thus, it becomes clear that the Newtonian analysis applies only 
to late times in the lifetime of our universe. For earlier times and larger scales one needs to 
employ general relativity. Curiously enough, it was not until 1957, long after Lifshitz's fully 
relativistic treatment, that Bonnor employed the Newtonian theory to study perturbations in a 
dust dominated FRW cosmology [^. In some ways the relativistic approach is simpler than the 
Newtonian, which requires considerable mathematical subtlety. For an extensive discussion of 
the Newtonian approach, we refer the reader to Peebles (1980) Q. Here, we will apply the idea 
of the Jeans instability to the simple case of an expanding self-gravitating non-relativistic fluid. 

2.1 The general fluid equations 

To begin with let us consider some fundamental ideas that apply to both Newtonian and rela- 
tivistic settings. We model the universe as a fluid, so that all the relevant quantities are described 
by smoothly varying functions of position. Cosmic strings, domain walls and other topological 
defects have no place in this picture. After the first bound structures, such as galaxies, form they 
are treated like particles along with the genuine particles that remain unbound. Note that such 
a fluid description of the universe applies only after smoothing over large comoving scales. An 
additional important concept is that of the "comoving observer" . Loosely speaking, a comoving 
observer follows the expansion of the universe including the effects of any inhomogeneities that 

^Throughout these lectures we will de dealing with cosmological models where the Hubble radius and the 
particle horizon, that is the maximum proper distance travelled by a typical photon, are effectively identical. We 
will therefore use the two concepts interchangeably. The reader is referred to Padmanabhan (1993) |^ for an 
illuminating discussion on the differences between these two important cosmological scales. 



3 



may be present. Adopting a Newtonian reference frame, specified by Cartesian space coordi- 
nates Tq, and a universal time t, we consider a fluid with density p and pressure p, moving with 
velocity in a gravitational potential Its evolution is governed by the standard Eulerian 
equations for a self gravitating medium 



dvn 



^ + dUpvn = 0, (1) 

„^ +v'^df3Va + -daP + da<^ = 0, (2) 

ot p 

d'^^-AirGp = 0. (3) 

where 9^ = d°'da is the associated Laplacian operator. Expressions (|^)-(|^ are respectively 
known as the continuity, the Euler and the Poisson equations. They describe mass conservation, 
momentum conservation and the Newtonian gravitational potential. The simplest solution to 
the above system is not applicable to cosmology, as it corresponds to a static matter distribution 
with p ,p = constant .0 For an expanding fluid it is convenient to adopt a "comoving" coordinate 
set (xa) instead of the "physical" (or proper) coordinates (ra) employed in The two 

frames are related via 

ra = axa , (4) 
where a = a{t) is the scale factor of the universe. The above immediately implies the relation 

Va = HVa + Ua , (5) 

between the physical velocity Va = dr^/dt and the "peculiar" velocity = dxQ,/dt, where 
H = a/a represents the Hubble parameter. In a comoving frame Eqs. read 

% + 2>Hp + -da{pu^) = 0, (6) 
ot a 

-T^Xa + + HUa + -U^dpUa + —daP + = , (7) 

dt^ ot a ap a 

d^^-AirGa'^p = 0. (8) 

One arrives at the above from (|l|)-@ on using the transformation laws {d/dt)phys = {d/dt)com — 
Hx°'da and (9Q-)phys = (l/a)(9a)com between physical and comoving derivatives. 

2.2 The unperturbed background 

The simplest non-static solution to the system @-® describes a smoothly expanding (i.e. Ua = 
0), homogeneous and isotropic fluid (i.e. po = po{t) , po = po{t)). In particular, the unperturbed 

^Greek indices indicate Newtonian spaces and take the values 1,2,3, while Latin characters run between and 
3 and denote spacetime quantities. 

■'in such a configuration the gravitational force vanishes (i.e. — from Eq. (^)). This clearly contradicts 
the Poisson equation, as it seems to suggest that the matter density vanishes as well. Following Jeans, we assume 
that Eq. (M) describes relations between perturbed quantities only (Jeans' swindle). 
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background universe is characterized by the system 

0, (9) 

0, (10) 
0, (11) 

with solutions 

po oc , Vq = Hva and <I>o = ^ttGo^pox'^ ., (12) 

where the expression for the gravitational potential follows from the isotropy assumption, namely 
from the fact that da = d/dxa = d/dx, and equation d^a/dt^ = —AnGap/S. Note that we no 
longer need to invoke Jeans' swindle, although the approach is still not entirely problem free, 
as $0 has a spatial dependence despite the assumption of spatial homogeneity. Also, on large 
scales Va > 1/ H \h, which substituted into Eq. (|l2|b) gives an expansion velocity greater 
than that of light. We remind the reader, however, that the Newtonian treatment applies to 
sub-horizon scales only. 



at'' a 
52$o-47rGaVo = 



2.3 The linear regime 

Consider perturbations about the aforementioned background solution 

p = po + 5p, p = po + 6p, t;" = t;^ + Jt;" , $ = «>o + 5$ . (13) 

where 5p, 5p, 5v" and 5^ are the perturbed first order variables with spatial as well as temporal 
dependence (i.e. 5p = 5p{t,Xa))- In the linear regime the perturbed quantities are much smaller 
than their zero order counterparts (i.e. 6p <^ po) During this period higher order terms, for 
example the product 6p6va, are negligible. This means that different perturbative modes evolve 
independently and therefore can be treated separately. Note that 6va = Ua is simply the peculiar 
velocity describing deviations from the smooth Hubble expansion. Also, the fluid pressure is 
related to the density via the equation of state of the medium. For simplicity, we will only 
consider "barotropic" fluids with p = p{p). 

Substituting (p^) into Eq. (P) and keeping up to first order terms only we obtain 

^-l + mdp + ^djv'' = (14) 

ot a 

on using the zero-order expression dp^/dt + ^iHpQ = (see Eq. (P)). The above describes the 
linear evolution of density fluctuations. In perturbation analysis, however, it is advantageous to 
employ dimensionless variables for first order quantities. Here, we will be using the dimensionless 
"density contrast" 6 = Sp/pQ. Throughout the linear regime 5p <^ po meaning that 5 <C 1. 
Introducing 6 we may recast Eq. (|T^ ) into 

^ + -daSv" = 0, (15) 
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given that 85 /dt = {1/ po)d5p/dt + {3H/ po)Sp from Eq. 1^. Also, substituting (|l3|) into Eq. (g) 
and linearizing we obtain the propagation equation for the velocity perturbation|j 

'd!!^ + H5v^ + ^d^5 + -do,5^ = Q. (16) 
at apQ a 

In deriving the above we have used the zero-order part of Eq. (^) and the binomial law {l+5)~^ ~ 
1 — 5 (recall that (5 <C 1). Also, for a barotropic fluid p = p{p) and daSp = v^daSp, where 
Vg = dp /dp is the adiabatic sound speed. Finally, inserting ( [l3| ) into the Poisson equation and 
linearizing we arrive at the relation 

9^(5$ - 47rGaVo<5 = , (17) 



for the evolution of the perturbed gravitational potential. Results (p~5[) -(P?7[) determine the 
behavior of the density perturbation completely. When combined they lead to a second order 
differential equation that describes the linear evolution of the density contrast. In particular, 
taking the time derivative of ( |T5|) and employing Eqs. ([T^), (^) we find that 

to linear order. We have arrived at the above by assuming Newtonian gravity, which requires 
matter domination and zero cosmological constant. As long as we remain well within the horizon, 
however, Eq. (|18|) also describes ordinary matter perturbations in the presence of a smooth 
radiative background or a cosmological constant. This is so because general relativity always 
reduces to Newtonian gravity near a free-falling observer moving at non-relativistic speeds. The 
only difference (see Peebles (1980) Q) is that the Poisson equation has been replaced by Eq. (|4l|). 
This changes the smooth gravitational potential but not its perturbation, provided the latter 
is due to matter alone. As a result, Eq. (|l^) is still valid if p is the simply density of the 
non-relativistic matter and 5 the associated density contrast. 

2.4 The Jeans Length 



Equation (18) is a wave- like equation with two extra terms in the right-hand side; one due to 
the expansion of the universe and the other due to gravity. It is therefore, natural to seek plane 
wave solutions of the form 



(k) 



where 5^y^ = (5(k)(t), k^ is the comoving wavevector and k = \/kak" is the associated comoving 
wavenumber.|^ Fourier decomposing Eq. ( |l8| ) and using the relations daS = ik^f^, d^5 = —\?5 
we obtain 



^Velocity perturbations decompose into rotational modes 5v'^, where daSvf^ = 0, and irrotational modes Sv^ 
with curWuy" = 0. However, it is only the divergence of Sv" (i.e. daSv") which contributes to the evolution of the 
density contrast (see Eq. (p^). Thus, rotational modes do not couple to the longitudinal density perturbations 
we will be addressing in these lectures. 

®The comoving wavelength of the perturbative mode is given by ^ = 27r/k, while A = is the physical (proper) 
wavelength. Also, the physical wavenumber is n = 2tt/X with n = k/a. Finally, we note that k^a;" — nar°'. 



6 



which determines the evolution of the k-th perturbative mode. The first term in the right-hand 
side of the above is due to the expansion and always suppresses the growth of Sq^^ . The second 
reflects the conflict between pressure support and gravity. When 'iirGpQ » v^k^ ja} gravity 
dominates. On the other hand, pressure support wins if v^k"^ /a? 3> A-kGpq. The threshold 
AnGpQ = vlk"^ /a? defines the length scale 

h=v.^. (21) 

The physical scale A j , known as the "Jeans length" , constitutes a characteristic feature of the 
perturbation. It separates the gravitationally stable modes from the unstable ones. Fluctua- 
tions on scales well beyond Aj grow via gravitational instability, while modes with A ^ Aj are 
stabilized by pressure]^ The Jeans length corresponds to the "Jeans mass" , defined as the mass 
contained within a sphere of radius Aj/2 

Mj = |vrp(^)', (22) 
where p is the density of the perturbed component. 



2.5 Multi-component fluids 



Thus far we have only considered a single-component fluid. When dealing with a multi- 
component medium (e.g. baryons, photons, neutrinos or other exotic particles), perturbations 
in the non-relativistic component evolve according to 



d^ 

dt2 



-2H 



d^ 



+ 



-Si 



(23) 



where the index i refers to the component in question. The sum is over all species and = 
Pi/po provides a measure of each component's contribution to the total background density 
po = YliiPi- Note that any smoothly distributed matter does not contribute to the right- 
hand side of the above. However, the unperturbed density of this component contributes to 
the background expansion. To first approximation, H is determined by the component that 
dominates gravitationally, while Vs is the velocity dispersion of the perturbed species which 
provide the pressure support.^ Equation ( |23| ) applies to the k-th perturbative mode, although 
we have omitted the wavenumber (and the tilde) for simplicity. 



2.6 Solutions 



We will now look for solutions to Eqs. (^), (23) in the following four different situations: 
(i) Perturbations in the dominant non-relativistic component (baryonic or not) for t > ieq- 

^One also arrives at the Jeans length following a simple qualitative argument. The timescale for gravitational 
collapse is fgrav ~ , whereas the timescale for the pressure forces to respond is tprcss ~ A/us. When 

tgrav <C fprcss, that is whcu A S> Ws(l/G'p)^^^ ~ Aj, pressure gradients do not have the time to respond and restore 
hydrostatic equilibrium. In the opposite case, namely for tprcss <Si tgrav, pressure totally overwhelms gravity. 

*The pressure of a baryonic gas is is the result of particle collisions. For "dark matter" , collisions are negligible 
and pressure support comes from the readjustment of the orbits of the colisionless species. In both cases it is the 
velocity dispersion of the perturbed component which determines the associated Jeans length. 
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(ii) Fluctuations in the non-baryonic matter for t < t^q- (iii) Baryonic perturbations in the 
presence of a dominant colhsionless species, (iv) Perturbations in the matter component during 
a late-time curvature dominated regime. 



2.6.1 Perturbed Einstein - de Sitter universe 

Consider a dust dominated (i.e. p = = v^) FRW cosmology with flat spatial sections (i.e. 17 = 
1). This model, also known as the Einstein-de Sitter universe, is though to provide a good 
description our universe after recombination. To zero order a oc t"^^^, H = 2/3t and po = 
l/GnGt'^. Perturbing this background, we look at scales well bellow the Hubble radius where 
the Newtonian treatment is applicable. On using definition (^l|) and the relation A = 27ra/k, 
Eq. (M) reads 



, 4,^ _ 2 
dt2 +3^dt 5 



6 = 0. (24) 



Note that form now on we will drop the tilde (~) and the wavenumber (k) for convenience. For 
modes well within the horizon but still much larger than the Jeans length (i.e. when Aj <C A ^ 
Xh), we find 

5 = Cit^l^ + C2t-^ , (25) 

for the evolution of the density contrast. As expected, there are two solutions: one growing and 
one decaying. Any given perturbation is expressed as a linear combination of the two modes. 
At late times, however, only the growing mode is important]^ Therefore, after matter radiation 
equality perturbations in the non-relativistic component grow proportionally to the scale factor 
(recall that a oc ^1"^ for dust). Note that baryonic perturbations cannot grow until matter has 
decoupled from radiation at recombination (we always assume that teq < ^rec)- Dark matter 
particles, on the other hand, are already colhsionless and fluctuations in their density can grow 
immediately after equipartition. After recombination, perturbations in the baryons also grow 
proportionally to the scale factor. 

On scales well bellow the Jeans length (i.e. for A ^ Aj), Eq. (p3) admits the solution 



b = Ct-l/6e±V2/3{Aj/A)lnt ^ (26) 

which describes a damped oscillation. Thus, small-scale perturbations in the non-relativistic 
matter are suppressed by pressure. Note that Aj/A oc i^/H, where u = Vs/X is the oscillation 
frequency. For Aj ^ A there are many oscillations within an expansion timescale (adiabati- 
cally slow expansion). Also, before recombination baryons and photons are tightly coupled and 
(fg)b cx Tb ~ oc a~^, implying that Aj/A = constant. After decoupling {v'^)h oc Tb oc a"^. 



2.6.2 Mixture of radiation and dark matter 

Consider the radiation dominated regime when a oc i^/^ and H = l/2t. On scales much smaller 
than the Hubble radius we can employ the Newtonian theory to study perturbations in the non- 
relativistic matter. The Newtonian equations are still applicable provided that the expansion is 

^Solution ( p5| ) demonstrates the difference between the Jeans instabiUty in the static regime Q (e.g. within 
a galaxy) and in the expanding universe. The expansion slows the exponential growth of the static environment 
down to a power law. 
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determined by the dominant radiative component. Applying Eq. ( p3D to a mixture of radiation 
and collisionless particles, with Vs = 0, we have 

d%M ^ 2H^^ = AttGpo {e^5^ + eoM^DM) ^ AirGpyd^ , (27) 

since pDM Pj — Po (i-e- ^dm 1 whereas ~ 1). Given that prior to equipartition H = l/2t 
and the small-scale photon distribution is smooth (i.e. {6y) ~ see Sec. 3.4.2), the above reduces 
to 

and admits the solution 

Sdm =Ci+ Cslnt . (29) 

Thus, in the radiation epoch small scale perturbations in the collisionless component experi- 
ence a very slow logarithmic growth, even when A > Aj.^ The stagnation or freezing- in of 
matter perturbations prior to equilibrium is generic to models with a period of rapid expansion 
dominated by relativistic particles and is sometimes referred to as the Meszaros effect §. 

2.6.3 Mixture of dark matter and baryons 

During the period from equilibrium to recombination perturbations in the dark component 
grow by a factor of Orec/oeq = Teq/Trec — 21Uh^ . At the same time baryonic fluctuations 
do not experience any growth because of the tight coupling between photons and baryons. 
After decoupling, perturbations in the ordinary matter will also start growing driven by the 
gravitational potential of the collisionless species. To be precise, consider the post-recombination 
universe, with a oc 

t2/3 and H = 2/3t, dominated by non-baryonic dark matter. Following 
Sec. 2.6.1, perturbations in the collisionless component grow as (5dm = Ca, where C is a constant. 



For baryonic fluctuations on scales larger that Aj, Eq. (23) gives 

^^-^ + = AnGpo (eoMf^DM + eb^^b) - 47rGpDM<5DM , (30) 

dt^ at 

since — for both species and pb Pdm — Po ■ Introducing the scale factor as the independent 
variable we may recast the above as 

where we have also used the relation a cx t^/^. The initial conditions at recombination are 5b = 0, 
because of the tight coupling between the baryons and the smoothly distributed photons, and 
^DM 7^ given that the dark matter particles are already collisionless. The solution 



shows that 6^ — > 6dm as a 3> a^ec- In other words, baryonic fluctuations quickly catch up with 
perturbations in the dark matter component after decoupling. Alternatively, one might say that 
the baryons fall into the "potential wells" created by the collisionless species. 



A more careful treatment, keeping the right-hand side in Eq. ( |27] ) gives 5dm cx: l-|-3a/2oeq (see Padmanabhan 
(1993); Coles & Lucchin (1995) Q). Hence, 5dm ~ constant for a a^q, namely during the radiation era. As 
the dust era progresses a 3> acq and 5dm oc a in agreement with Sec. 2.6.1. 
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2.6.4 The curvature dominated regime 



So far, and also for the rest of these notes, we have assumed a high density background universe 
with cfose to unity. If, instead, is small the universe could become curvature dominated 
at late times. Soon after the transition into the curvature regime, the rapid expansion of the 
universe prevents the inhomogeneities from growing. Let us take a closer look at this issue. 
When curvature dominates <C 1, a (x t and H = 1/t (see Kolb & Turner (1990) H). For 



collisionless matter with Vs — Eq. (20) gives 



*=g + 2*i^ = in.. (33) 



since AirGpo = 3^2/2^^. Moreover, given that 0^1 the right-hand side of the above is effectively 
second order, and Eq. (|33|) reduces to 

The solution 

5 = Ci+C2t"S (35) 
verifies that perturbations cease growing when curvature dominates. 



2.7 Summary 

The Newtonian treatment suffices on sub-horizon scales and as long as we we deal with fluc- 
tuations in the non-relativistic component. During the radiation and the curvature dominated 
eras, perturbations do not grow. This suppression occurs because, in both epochs, the expansion 
is too rapid for the perturbations to experience any growth. After equipartition, fluctuations 
with wavelengths larger than the Jeans length grow as 5 oc a oc t^/^, whereas those bellow Aj 
oscillate like acoustic waves. Note that baryonic perturbations do not grow until recombination 
due to the tight coupling between baryons and photons. Dark matter fluctuations, on the other 
hand, start growing immediately after matter-radiation equality. As soon as the baryons have 
decoupled perturbations in their density distribution will be driven by the gravitational poten- 
tial of the collisionless species. Thus, shortly after recombination, baryonic fluctuations grow 
rapidly and soon equalize with those in the dark matter. Subsequently, perturbations in both 
components grow proportionally to the scale factor. 



3 Linear relativistic perturbations 

Although the Newtonian analysis provides valuable insight into the behavior of inhomogeneities, 
it has serious shortcomings. The proper wavelength of any perturbative mode will be bigger 
than the horizon at sufficiently early times. On such scales general relativistic effects become 
important and gauge ambiguities need to be settled. Moreover, one cannot use the Newtonian 
theory to study perturbations in the relativistic component. Given that all the astrophysically 
relevant scales were outside the horizon early on (i.e. for z > 10^), it becomes clear that a general 
relativistic treatment of cosmological density perturbations is imperative. General relativity was 
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first applied to cosmological perturbations in a seminal paper by Lifshitz in 1946 @. The Lifshitz 
approach relies upon selecting a "gauge" , finding the solutions in that gauge and then identifying 
the "gauge modes" . In the relativistic treatment one perturbs both the spacetime metric and the 
energy-momentum tensor of the matter sources. In other words, one assumes that gab = gab~^^9ab 
and Tab = T°b ~^ ^"^ab: where the zero index denotes the background variables. For "small" 6gab 
and 5Tab, one can perturb and then linearize the Einstein Field Equations (EFE). Here, we 
will only outline the key steps of the analysis and state the main results, while we refer the 
reader to ^ for the details. Note that the complete set of the relativistic equations reveals 
three types of perturbations, namely tensor, vector and scalar modes. Tensor perturbations 
correspond to the traceless, transverse part of 6gab- They describe gravitational waves and have 
no Newtonian analogue. Vector and scalar perturbations, on the other hand, have Newtonian 
counterparts. Vector modes correspond to rotational perturbations of the velocity field, while 
scalar modes are associated with longitudinal density fluctuations. In these lectures we will only 
consider the latter type of perturbations. 

3.1 The gauge problem 

It has long been known (see Lifshitz (1946) Q) that the study of cosmological perturbations is 
plagued by the so called gauge problem, which stems from the fact that in perturbation theory 
we actually deal with two spacetimes.]^ The physical spacetime, W, corresponding to the real 
universe and a fictitious one, W, which defines the unperturbed background. The latter is usually 
represented by the homogeneous and isotropic FRW models. To proceed we need to establish 
an one-to-one correspondence 99 : W i— > W, namely a gauge, between the two spacetimes. When 
a coordinate system is introduced in W the gauge caries it to W and vice versa, thus defining 
a background spacetime into the real universe. In practice what a gauge does is define the 
"slicing" and "threading" of the spacetime into spacelike hypersurfaces and timelike worldlines 
respectively. Any change in : W 1— > W, keeping the background coordinates fixed, is known 
as a "gauge transformation". The latter induces a coordinate transformation in the physical 
spacetime but also changes the event in W that is associated with a given event in W. Therefore, 
gauge transformations should be distinguished from coordinate transformations which simply 
relabel events.]^ The problems stem from our freedom to make gauge transformations. By 
definition, a perturbation in a given quantity is the difference between its value at some event in 
the real spacetime and its value at the corresponding event in the background. This means that 
even if a quantity behaves as a scalar under coordinate changes, its perturbation will not be 
invariant under gauge transformations, provided that the quantity in question is non-zero and 
time dependent in the background. As a result, one might end up with spurious gauge modes 
in the solutions, which have no physical meaning whatsoever. Typical example are density 
perturbations. Given that the density is a time dependent scalar in the background, density 
perturbations are not invariant under gauge transformations that change the correspondence 
between the hypersurfaces of simultaneity in W and W. Note that on scales bellow the horizon 
the hypersurfaces of constant time are physically unambiguous and the gauge choice is well 

^^For a extensive discussion on the gauge problem in cosmology the interested reader is referred to articles by 
Sachs (1966) @, Bardeen (1980) (ll), Ellis & Bruni (1989) ||l|, Stewart (1990), Bruni et al (1992) ||l^ 

^^In practice, we often talk about coordinate changes in W, since a frame in W corresponds, via (p : W 1-^ W, 
to one already established in W. Therefore, the coordinate choice in W also determines the gauge between W 
and W. In this respect, gauge transformations can be represented as coordinate changes in W. 
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motivated. On super-horizon scales, however, one is free to choose between gauges that give 
entirely different results for the time dependence of the density perturbation. It becomes clear, 
therefore, that Newtonian perturbations are not plagued by gauge ambiguities but early universe 
studies are. 

There are two ways of dealing with the gauge problem. The first, which will be our approach, 
is to choose a particular gauge and compute everything there. If the gauge choice is well 
motivated, the perturbed variables will be easy to interpret. However, the task of selecting the 
best gauge for any given situation, something known as the "fitting problem" in cosmology, 
is not always trivial. The second approach is to describe perturbations using gauge-invariant 
variables. In a very influential paper, Bardeen (1980) Q introduced a fully gauge- invariant 
approach following earlier work by Gerlach & Sengupta (1978) (see also Kodama &: Sasaki 
(1984) for an extensive review). Bardeen's approach, however, is of considerable complexity as 
it determines a set of gauge-invariant quantities that are related to density perturbations but are 
not perturbations themselves. Building on earlier work by Hawking (1966), Stewart & Walker 
(1974) and Olson (1976), Ellis & Bruni (1989) |12| formulated a fully covariant gauge-invariant 
treatment of cosmological perturbations. Their approach, which is of high mathematical elegance 
and physical transparency, has the additional advantage of starting from the fully non-linear 
equations before linearizing them about a chosen background. 



3.2 The relativistic equations 

We proceed by adopting what is usually known as the covariant Lagrangian approach to cos- 
mology, which is also a direct extension of the Newtonian treatment. The basic philosophy is 
to employ locally defined quantities and derive their evolution equations along the worldlines of 
the comoving observers. We start by assuming a congruence of timelike worldlines tangent to 
the 4- velocity vector Ua- The later determines the motion of a fundamental observer comoving 
with the fluid and is normalized so that uq = 1 = u'^ and = 0, with a = 1,2,3. Note that 
throughout these lectures we have set c = 1. If gab is the spacetime metric, with signature 

(H ), then hab = gab — UaU^ is the metric of the three-dimensional spaces orthogonal to the 

observer's motion (note that habU^ = 0), which define the observer's instantaneous rest space. 
Also, if V(j is the covariant derivative relative to gab, then = ha^^b is the covariant derivative 
operator on these hypersurfaces provided that there is no vorticity. Also, = h^^VaVb is the 
associated Laplacian. Finally, an overdot indicates differentiation along Ua, namely derivatives 
with respect to proper time r. For example, Ua = u^Vb^a is the 4-acceleration. The above 
procedure determines our gauge choice. 

The basis of the relativistic analysis is the Einstein Field Equations (FEE) describing the 
interaction between matter and spacetime geometryf^ 

Gab = Rab — \R9ab = ^^afc " ^Qab , (36) 

where Gab is the Einstein tensor, Rab is the Ricci tensor, R = Ra"" is the Ricci scalar. Tab is the 
total energy-momentum tensor of the mater fields, A is the cosmological constant and n = SirG. 
The Einstein tensor has the extremely important property of having an identically vanishing 
divergence, that is V'^Gab = 0. When applied to Eq. (|^), the latter yields the conservation law 

y'Tab = 0. (37) 

^''For an updated and extensive discussion on relativistic cosmological models the reader is referred to Ellis & 
van Elst (1999) (l3|. 
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For a perfect fluid the stress-energy tensor takes the simple form 

Tab = PUaUb + Phab , (38) 

where p and p are respectively the energy density and pressure of the fluid. Here, similarly to the 
Newtonian analysis, we assume a barotropic fluid with p = p{p). Note that the assumption of an 
isotropic pressure is not actually valid before matter-radiation equality owing to particle diffusion 
and free streaming effects (see Sees. 4.5.1 and 5.4). The relativistic analogues of the continuity 
and Euler equations are obtained from the timelike and spacelike parts of the conservation law 
(|37|). In particular, substituting (|^) into Eq. (|37|) and then projecting along the observer's 
motion (i.e. contracting with Ua) we obtain the energy density conservation law 

p + ?,H{p + p) = {), (39) 

where 3if = 'V"'Ua = D^Ua. On the other hand, by projecting orthogonal to Ua we arrive at the 
momentum density conservation equation 

{p + p)Ua + 'DaP = 0. (40) 

Equations (^), (^0|) are supplemented by 

RabU^u^ = ^k{p + 3p) . (41) 

The above relates the spaccetime geometry to the matter sources and is the relativistic analogue 
of the Poisson equation. Note that throughout these notes the cosmological constant has been 
set to zero. 



3.3 The linear regime 

Before perturbing Eqs. (p9|), (pO|), we need to make one additional step, which will fix our 
gauge completely. So far we have not assigned time labels to the comoving hyper surf aces. 
The proper time interval r is position dependent on these surfaces and does not provide an 
appropriate label for coordinate time. Assuming that t is a valid ordering label, one can show 
(see e.g. Padmanabhan (1993), Liddle & Lyth (2000) §) that 

1 = 1-^. (42) 

at P + P 

On using the above, the perturbed continuity equation gives 

{SpY = -3{po + po)SH - 3Ho6p (43) 

to first order, where the dash indicates differentiation with respect to t. Defining 6 = 5p/pQ the 
above is recast into 

6' - 3wHo6 + 3(1 + w)6H = , (44) 



where w = po/po determines the equation of state of the medium. Also, starting from Eq. (|4 
we have 



{6Hy + 2Ho6H + l7rGpo6 + —^DH = 0, (45) 

3(1 + w) 
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to first order, where 5H describes scalar deviations from the smooth background expansion. We 
obtain Eq. ( p5|) by taking the 4-divergence of Eq. (|40| ) and then employing the Ricci identity. 
Applied to the 4- velocity vector the latter reads V[aVfe]Uc = Rahcd.u'^: where Rahcd is the spacetime 
Riemman tensor {Rah = R'^ach is the associated Ricci tensor). To proceed further we need the 
following auxiliary formulae 

w' = 3Ho{l + w){w- v^) and H'q = -^{1 + w)H^ , (46) 

where the latter is commonly referred to as the Raychaudhuri equation. Taking the time deriva- 
tive of Eq. ( ^5| ) , substituting (^) , using the auxiliary relations ( ]4^ ) and keeping up to first order 
terms we arrive at 

1,2 2 

+ {2-6w + 3z;2) Ho~5[^^ _ 3 (i + _ 3^2 _ g^2^ ^2^^^^ ^ "^^(k) ' (^7) 

for the evolution of the k-th perturbative mode. In deriving the above we have also employed 
the Fourier decomposition 6 = ^i^\k)Q^^^ ^oi the density perturbation, with Da^(k) = 0. Here 
Q(k) are scalar harmonics, with Q'q^-^ = 0, which are solutions of the Laplace-Beltrami equation 

D'Q(k) = -^Q(k)- (48) 

The above immediately implies that = — (k^/a^)5, which explains the term in the right- 
hand side of Eq. ( p7|) . This equation may be thought of as the relativistic counterpart of the 
Newtonian formula (|20|) . One recovers the Newtonian limit by setting w = = v'^ in the left 
hand side of (^7|), and using the background Friedmann equation Hq = 8irGpo/3. 

Note that in many applications it helps to recast Eq. (47) in terms of the scale factor. 
Introducing a as the independent variable we find 



'l + 8w- 3u;2 - 6t;,2 



^(k)=0, (49) 



where we have employed the transformation laws d/dt = aH d/da and d^/dt^ = a^H^ d^/da^ 
[{l + 3w)aH^/2] d/da. 



3.4 Solutions 

We will now seek solutions to the relativistic perturbation equations to supplement the New- 
tonian results of the previous section. The cases to be considered are: (i) Super-horizon sized 
perturbations in the dominant non-relativistic component after equilibrium, (ii) Fluctuations 
in the relativistic matter before matter-radiation equality both inside and outside the Hubble 
radius. 



3.4.1 Perturbed Einstein - de Sitter universe 

The Newtonian analysis is valid for modes well within the Hubble radius. On scales beyond Xh, 
however, one needs to engage general relativity even when dealing with non-relativistic matter. 



For pressureless dust w = = and Eq. (49) reads 



od^5 o d<5 o / 2 k2?;2 . 

+ i-i-o7^M = o, (50) 



da^"^ da H "^a^H'o 
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where again we have dropped the tilde and the wavenumber for simphcity. For modes lying 
beyond the Hubble radius A Xh and ]^v1 / o^Hq <^ 1. On these scales the above reduces to 



da"= ^ da 



with 

6 = Cia + Caa^^/^ . (52) 
Thus, after recombination large-scale perturbations in the non-relativistic component grow as 

3.4.2 The radiation dominated era 

Before equipartition radiation dominates the energy density of the universe and w = 1/3 = v^. 
During this period Eq. ( ]49| ) gives 

On large scales, when A » Xh and k^/o^i^g <^ 1, the above reduces to 

^2^_2A =0, (54) 

with a power law solution of the form 

6.y = Cio^ + 620,^^ , (55) 

where Ci, C2 are constants. Hence, before matter-radiation equality large-scale perturbations 
in the radiative fluid grow as 6^ oc a^/^. Note that Eq. (|5^ ) also governs the evolution of the 
non-relativistic component (baryonic or not), since it does not incorporate any pressure effects. 
Therefore, solution ( |55[) also applies to baryons and collisionless matter. 
On sub-horizon scales, with A <C A^ and k'^/a^HQ » 1 Eq. ( |53| ) becomes 

and admits the oscillatory solution 

6^ = Ce'^"/^, (57) 

where C = constant. Thus, in the radiation era small-scale fluctuations in the relativistic com- 
ponent oscillate like sound waves. Given that Xh/X ^ 1, the oscillation frequency is very high. 
As a result, (5^) ~ on scales well bellow the Hubble radius. In other words, the radiative fluid 
is expected to have a smooth distribution on small scales (see Sec. 2.6.2). 

Note that in the radiation era the transition from growing to oscillatory modes occurs at 
A ~ Xh, which implies that before equipartition the role of the Jeans length is played by the 
Hubble radius. 
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3.5 Summary 



General relativity is necessary on scales outside the horizon and also when studying perturbations 
in the relativistic component. Prior to equilibrium, fluctuations in the radiative fluid grow as 
djona'^cxit for wavelengths larger the Hubble radius. Super-horizon sized matter perturbations 
also grow at the same rate. On small scales, however, fluctuations in the photon density oscillate 
rapidly. As a result, the small-scale distribution of the radiative fluid, and of the tightly coupled 
baryons, remains smooth. After equipartition perturbations in the non-relativistic matter grow 
proportionally to the scale factor. 

4 Baryonic structure formation 

We will now consider cosmological models where baryons are the dominant form of matter. It 
should be made clear at the outset, however, that purely baryonic models cannot successfully 
explain the origin of the observed structure. Nevertheless, it is important to look at the details 
of these scenarios. After all, whatever the dominant form of matter, baryons do exist in the 
universe and it is crucial to study their behavior. The key issue is understanding the interaction 
between baryonic matter and radiation during the plasma epoch. The simplest way of doing so 
is by looking at models containing these two components only. 

4.1 Adiabatic and isothermal perturbations 

Before recombination the universe was a mixture of ionized matter and radiation interacting 
via Thomson scattering.]^ When dealing with the pre-recombination plasma we distinguish 
between two types of perturbations, namely between "adiabatic" and "isothermal" disturbances 
(Zeldovich (1967) |l^). The former include fluctuations in both the radiation and the matter 
component (i.e. dhj^-y 0)) whereas in the latter only the matter component is perturbed 
(i.e. 5b 7^ but 6-y = 0). Before recombination, a generic perturbation can be decomposed into 
a superposition of independently propagating adiabatic and isothermal modes. After matter 
and radiation have decoupled, however, perturbations evolve in the same way regardless of their 
original nature. Because there is no interaction between baryons and photons and radiation field 
is dynamically negligible, the universe behaves as a single-fluid dust model. 

4.1.1 Adiabatic perturbations 

Adiabatic (or isentropic) modes contain fluctuations both in the matter and the radiation com- 
ponents, while keeping the entropy per baryon conserved. Note that if 5 = 4py/3T^ is the 
photon entropy, then S = S/kBiih is the entropy per baryon, where T^ is the photon tempera- 
ture, /cb is the Boltzmann constant and Ub is the baryon number density. Given that oc T^ 

and pb Of Ub, the entropy per baryon satisfies the relation S oc pS^^^/pb- Consequently one arrives 
at 




(58) 
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For simplicity we will neglect the presence of Helium nuclei and the role of the neutrinos. 
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which is known as the condition for adiabaticity. A set of density contrasts satisfying the 
above requirement consists an adiabatic perturbation. The later are naturahy generated in the 
simplest inflationary models through the vacuum fluctuation of the inflaton field (see Liddle & 
Lyth (2000) @). 



4.1.2 Isothermal and isocurvature perturbations 

Isothermal modes include only matter fluctuations, while the radiation field is assumed to be 
uniformly distributed. This means that the radiation temperature is also uniform (recall that 
p-^ oc T^), which explains the name isothermal. This type of perturbation is closely related 
to the "isocurvature" modes, where (5b / but 5p = Q [p = p\^ + p-y). This implies that 
the geometry of the 3-dimensional spatial hypersurfaces remains unaffected, hence the name 
isocurvature. Note that isocurvature perturbations can still affect, via the pressure, the 4- 
dimensional spacetime geometry. For such modes we have 

5p = ^ Pb^b + p^5y = Q ^ -L = -P^. (59) 

Oh P-y 

Clearly, when p^ ^ pb, as it happens in the radiation era, one finds that <C (5b, which explains 
why isocurvature modes are very often referred to as isothermal. Unlike adiabatic disturbances, 
isocurvature perturbations are usually absent from the simplest models of inflation. They can 
still be produced, however, in multi component inflationary scenarios by the vacuum fluctuation 
of a field other than the inflaton (see Liddle &: Lyth (2000) [Q]). As stated earlier, the distinc- 
tion between adiabatic and isothermal fluctuations is meaningful only prior to recombination, 
when matter and radiation are tightly coupled. After decoupling, perturbations in the matter 
component evolve as if they were effectively isothermal. 



4.2 Evolution of the sound speed 

The different nature of the adiabatic from the isothermal perturbations means that each type of 
disturbance has its own discrete signatures. The key issue is the evolution of the sound speed, 
since it determines the scale of gravitational instability. Here, we consider the evolution of the 
sound speed for the adiabatic and the isothermal modes during the plasma era. 



4.2.1 The adiabatic sound speed 

In a mixture of radiation and matter the total density and pressure are p = p-y + ph and 

p 1^ Py = Py/2> respectively (recall that p^ = Py/2> and pb — 0). Hence, the adiabatic sound 
speed is given by 



dp) ^\ Ap 



where we have used the adiabatic condition dp]^/dpy = Spy^/Ap^ (see Eq. (|5^)). In the radiation 

fa) r- 

era p-y » ph ensuring that i)s ~ l/v3. In the interval between equipartition and decoupling, 
when pb ^ P7! Eq. ( |60[ ) gives vi^^ ~ Y^4p^/3pb ^ a~^l'^ . In particular, employing the relations 
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[(1 + z)/{l + Zeq)]^ and Ph = Pcq[(l + + -^cq)]^, we find (see Coles & Lucchin 

(1995) i) 

vf^ - 2 X 10'^ 1^ ^ ' cm/sec. (61) 

Note that throughout these notes we assume that Zeq > z^^^, that is equipartition occurs prior 
to recombination. Given that 1 + Zgq ~ 4 x 10^(r2b^^) and 1 + z^-^c — 10^, our assumption holds 
only if n^h'^ > 0.02. 

4.2.2 The isothermal sound speed 

The sound speed associated with isothermal fluctuations is that of a monatomic gas 

where Tb is the matter temperature, is the proton mass and 7 = 5/3 for hydrogen. Be- 
fore decoupling photons and baryons are still tightly coupled and share the same temperature 
(i.e. Tb ~ T^). Equation ( |6^ ) then implies that vlf' oc T^^ oc a~"^/^. In particular, given that 
(T^)rec — 4 X 10^ and using the relation T^ = (T^)i.ec(l + + Zrec), we obtain (see Coles &: 
Lucchin (1995) |]) 

v'p = 5 X 10^ i^ l + z ) ^^/^^^ ' (63) 



for z < Zj-cc- After decoupling the isothermal sound speed obeys Eq. (63) as long as Tt, — T. 



■-7 

(i.e. for z > 300). Subsequently to that and until the time of reheating Tb ^ T^ and Vs oc 
T^^ oc a . 

As it was pointed out earlier, in the post-recombination we can no longer distinguish between 
adiabatic and isothermal perturbations. After decoupling the sound speed of matter perturba- 
tions coincides with the isothermal one. Thus, at recombination the sound speed of the adiabatic 
disturbances drops from vi^^ oc {p-y/ph)^^"^ to vi"^^ oc {ph/phy^^ (see Eqs. (|60|), (|62|)). Given that 
p^ ~ n^T^ and pb ~ iibTb with n^/ub ^ 10®, the reduction in vi^^ is very large (recall that 
Tb — T^ through recombination). Following such a drop in the adiabatic sound speed at decou- 
pling suggests, one anticipates a similarly dramatic decrease in the Jeans length at around the 
same time. 

4.3 Evolution of the Jeans length and the Jeans mass 
4.3.1 Adiabatic modes 

The difference in the sound speed between adiabatic and isothermal fluctuations means that 
their respective Jeans length and Jeans mass also differ. This in turn implies that adiabatic 
and isothermal modes become gravitationally unstable at different scales. In particular, using 
definition (|2l] ) and the results of Sec. 4.2, we find that before recombination the adiabatic Jeans 
length evolves as 

^(a) ^ I ^ > ^eq , ^g^^ 
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In deriving the above we have taken into account that vi'^^ ~ 1^1/3 and p ^ (x a ^ for 

z > Zeq, while vi^^ oc a~^/^ and p ~ oc when Zcq> z > Zj-^c- On using definition (^), the 
above result translates into 

[ constant Zcq > Z > Z^ec , 

with (Aff ^)i.cc =^ 3 X 10i5(17/i2)-2 M0 (see Coles & Lucchin (1995) Consequently, in the 

adiabatic scenario, the first scales to become gravitationally unstable and collapse soon after 
decoupling have the size of a supercluster of galaxies. 



4.3.2 Isothermal modes 

Throughout the plasma era Tb ^ and the isothermal sound speed evolves as Vg^ oc a^^/^, 
which substituted into definition ( |2l| ) gives 

A«oc|"'^' (66) 

The corresponding Jeans mass is governed by 

Mf ) oc I ^ > ^eq , .g7) 

[ constant Zeq > Z > Z^ec ■ 

with (MP)rec ~ 5 X 10^(17/i2)-i/2 Mq (sec Coles & Lucchin (1995) @). In the isothermal 
models the first sizes to collapse are of the order of a globular star cluster. According to results 
(pSD, (|67[), the Jeans mass, of both adiabatic and isothermal perturbations, increases during the 
radiation era but it remains constant in the interval between equipartition and recombination. 
So, Mj takes its maximum possible value in models with Zgq = -z^rec- Note the huge drop, of the 
order of 10^^ Mq, in M^""^ around decoupling. This is the result of a similarly large drop of the 
adiabatic sound speed at the same time (see previous section). 

After recombination the Jeans length and the Jeans mass of matter perturbations are taken 
to coincide with Aj"* and m]^^ respectively. The latter evolves as Mj oc a~'^/^, since Vs'^ oc a^^ 

for Z < Zrpc- 



4.4 Evolution of the Hubble mass 

An additional important scale for structure formation is that of the Hubble mass Mjj defined 
as the total amount of mass contained within a sphere of radius Ah/2, 

MH = l7rp(^^y, (68) 

where p is the density of the perturbed component and Xh is the Hubble radius. In these lectures 
we will only consider the baryonic Hubble mass (i.e. p = p\y always). Also, recall that we will only 
be dealing with models where the Hubble radius is effectively identical to the particle horizon. 

^^For consistency the values for the baryonic Jeans mass given in this section have all been quoted from Coles 
& Lucchin j^, with the assumption that f2 = S7b = 1. 
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In this respect, Xh and Mh define the scale over which the different parts of a perturbation 
are in causal contact. Note that a mass scale M is said to be entering the Hubble radius when 
M = Mfj. Given that Xh cetera? before equilibrium and Xh oc t oc a^/"^ for t > tgq we obtain 
the evolution law 

MhocI^L (69) 

It should be emphasized that before matter-radiation equality the baryonic Jeans mass is of the 
same order with the Hubble mass. Indeed, when radiation dominates p ~ p^, which means that 
Aj ~ (Gp^)~^/^ and Xh ~ {Gp^)~^^'^. Given that and using definitions ( p^ ) and ( p^ ) one can 
easily verify that Mj ~ Mh throughout the radiation epoch. 

4.5 Dissipative effects 

To this point we have treated the cosmic medium purely gravitationally, as a perfect fluid, 
ignoring any dissipative effects. In the process we have established two key physical scales, the 
Jeans mass and the Hubble mass, which play an important role in all structure formation models. 
We now turn our attention to other physical processes that can modify the purely gravitational 
evolution of perturbations. In baryonic models the most important physical phenomenon is 
the interaction between baryons and photons in the pre-recombination era, and the consequent 
dissipation due to viscosity and heat conduction. 

4.5.1 Collisional damping of adiabatic perturbations 

Adiabatic perturbations in the photon-baryon plasma suffer from collisional damping around the 
time of recombination because the perfect fluid approximation breaks down. As we approach 
decoupling, the photon mean free path increases and photons can diffuse from the overdense into 
the underdense regions, thereby smoothing out any inhomogeneities in the primordial plasma. 
The effect is known as collisional dissipation or "Silk damping" (Silk (1967) |]l5|). A detailed 
treatment requires integrating the Boltzmann equation through recombination. Here we will 
only obtain an estimate of the effect. To begin with, we consider the physical (proper) distance 
associated with the photon mean free path 

= ~ lO^^a^X-i i^hh^y^ cm , (70) 

AgneCTT 

where Xg is the electron ionization factor, Ug oc is the number density of the free electrons 
and (Tt is the cross section for Thomson scattering. Clearly, all baryonic perturbations with 
wavelengths smaller than £^ will be smoothed out by photon free streaming. The perfect fluid 
assumption breaks down completely when X <^ Damping, however, occurs on scales much 
larger than as the photons slowly diffuse from the overdense into the underdense regions, 
dragging along the still tightly coupled baryons. Within a time interval At a photon suffers 
N = At/i^ collisions and performs a random walk with mean square coordinate (i.e. comoving) 
displacement given by 

{Axf= NfM'= ^At, (71) 
\ a J 
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where {^^/a) is the coordinate distance between successive cohisions. Integrating the above up 
to recombination time we obtain the total coordinate distance travelled by a typical photon 



/■tree P ( \ /"tree 

_ 3 (■^7)rec^rec /^r)N 
~ 5 „2 ' V'^j 

"rec 

on using the relations = (^7)rcc(o/arec)^ and a = a^ccit/'tvecf'^^ for the photon mean free 
path and the scale factor of the universe respectively. The latter is a reasonable approximation 
given that earlier on the coupling between the photons and the electron was too tight for the 
photons to move at all. The physical scale associated with the above result is (see Kolb & Turner 
(1990) i) 



= axs = ^/|{l^UctZ ^ 3.5 {Qh^) Mpc , (73) 
and the associated mass scale, which is known as the "Silk mass", is given byf^ 

Ms = |m (I)' 6.2 X 10^2 [nh^)-'/^ Mo , (74) 

for adiabatic baryonic perturbations. Note that if we had not included dissipation, the amplitude 
of an acoustic wave on a mass scale smaller than the Jeans mass would have remained constant 
during the radiation era and then decayed as in the interval between equilibrium and 

decoupling (see Eq. (p6|)). The dissipative process we considered above causes the amplitude 
of these waves to decrease at a rate that depends on the size of the perturbation. The final 
result is that fluctuations on scales bellow the Silk mass are completely obliterated by the time 
of recombination and no structure can form on these scales. Alternatively, one might say that 
adiabatic perturbations have very little power left on small scales. 

4.5.2 Freezing-in of isothermal perturbations 

Consider isothermal perturbations in the pre-recombination era on scales larger than the Jeans 
mass. As the matter particles try to move around, they encounter a viscous friction force due to 
Thomson scattering with the smoothly distributed background photons. This acts an effective 
drag force on the baryons causing the isothermal mode to freeze-in throughout the plasma era. 
Qualitatively, one can explain this effect on the basis of the following simple physical argument. 
Consider the viscous force, per unit mass, due to the aforementioned "radiation drag" 

~ 7^ ~ , (75) 

where A is the wavelength of the perturbation and is the timescale for Thomson scattering. 
On the other hand, the self-gravitating pull per unit mass of a baryonic perturbation is 

Fg ~ — ~ GpbX ~ ^ , (76) 



^*'Given that tree ^ {£"/)icc we have {£^)Tectiec ^ {£j)iec, which ensures that ^ (^7)rcc- Consequently, Silk 
damping is felt on scales much larger than the mean free path of a typical photon. 
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since F ~ and /Ob 1/Gt^. Before recombination t^^ <<^t, implying that Fq <C -Ft- As a 
result isothermal perturbations cannot grow until matter and radiation have decoupled. Note 
that the stagnation of isothermal modes prior to recombination due to the radiation drag is not 
related to the Meszaros effect discussed in Sec. 2.6.2. The latter is purely kinematical and does 
not involve any interactions between matter and radiation. 

4.6 Scenarios and problems 

Depending on the nature of the primeval perturbations one may consider two scenarios of bary- 
onic structure formation. Those where the original inhomogeneities are of the adiabatic type 
and those permeated by isothermal fluctuations. The adiabatic and isothermal scenarios were 
in direct competition throughout the 1970s. One aspect of the confrontation was that the adia- 
batic models were advocated by the Soviet school of astrophysicists led by Zeldovich in Moscow, 
whereas the isothermal scenarios were primarily an American affair promoted by Peebles and the 
Princeton group. At the end, generic shortcomings in both models meant that neither of these 
adversaries won the battle. So in the early 1980s, the purely baryonic models were sidelined by 
non-baryonic dark matter scenarios (see Sec. 5). 

4.6.1 Adiabatic Scenarios 

Typically, adiabatic perturbations with sizes larger than the maximum value of the Jeans mass, 
that is with M > (Mj'^^)eq ~ 3x 10^^(ri/i^)~^ Mq experience uninterrupted growth. In particular 
they grow as 5b oc ^ before matter-radiation equality and like (5b oc t^/^ after equipartition. 
Fluctuations on scales in the mass interval (Mj''^)eq > M > (Ms) rec grow as (5b cx: t while they 
are still outside the Hubble radius. After horizon entry and until recombination these modes 
oscillate like acoustic waves. The amplitude of the oscillation is constant before equilibrium but 
decreases as between equipartition and recombination. After decoupling the modes become 
unstable again and grow as (5b oc t^^^. Finally all perturbations on scales smaller than the value 
of the Silk mass at recombination, that is with M < (Ms)rec — 6.2 x 10-*^^ (fi/i^) Mq are 
eventually dissipated by photon diffusion. In short, only fluctuations on scales exceeding that 
of a galaxy cluster can survive the plasma epoch. Schematically, the evolution of an adiabatic 
mode with mass scale M(A) is depicted in Figure 1. We distinguish between three evolutionary 
stages A, B and C, depending on the size of the perturbation and on the time of horizon crossing. 

In adiabatic scenarios the smallest scales with any structure imprint upon them at the time 
of recombination are as large as a rich cluster of galaxies. Fluctuations on smaller scales have 
been completely obliterated by Silk damping. After decoupling perturbations grow steadily 
until their amplitude becomes of order unity or larger. At that point the linear theory breaks 
down and one needs to employ a different type of analysis. Qualitatively, what happens is that 
those huge perturbations undergo anisotropic collapse to form massive flattened objects known 
as "pancakes" or "caustics" (Zeldovich (1970) ^^). After pancake formation, non-linear gas 
dynamics inside the collapsed structure lead to shock generation and the subsequent cooling 
causes fragmentation into smaller structures. Thus, in the adiabatic scenario galaxies are born 
out of larger condensations in a "top-down" fashion. The expected large-scale structure pattern 
is one of large sheet-like filaments with enormous voids in between, a picture that seems to fit 
qualitatively very well with observations. 

The main difficulty with the adiabatic scenario is that it predicts angular fluctuations in the 
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Figure 1: Typical evolution of a perturbed scale M ~ M{X) in adiabatic scenarios. During stage A 
{t < icnt < icq) the mode grows as (5b c>c t. Throughout stage B (tent < t < t^q) the perturbation 
oscillates. Finally, in stage C (tree < t) the mode becomes unstable again and grows as oc t^^"^ Note 
that fluctuations with size smaller than Ms are wiped out by photon diffusion. 



CMB temperature in excess with observational limits. More specifically, in a 17 = 1 baryonic 
universe, one needs 6h — 10"'^ at recombination if structure were ever to form by today. In 
the adiabatic picture, however, matter inhomogeneities are accompanied by perturbations in 
the radiation field. This will inevitably lead to temperature fluctuations of order 6T /T ~ 5^ ~ 
10^^ at decoupling, which is in direct disagreement with observations. To make things worse, 
primordial nucleosynthesis restricts the baryonic contribution to the total density of the universe 
down to Jib — 0.1. Given that perturbations grow slower in open models than in flat ones, a 
greater temperature fluctuation at recombination is required if structure were to form by now. 

4.6.2 Isothermal scenarios 

In isothermal models structure formation proceeds very differently. Isothermal fluctuations 
remain frozen-in throughout the radiation era due to the radiation drag. After recombination 
modes larger than (Mj''*)i.ec — 5 x 10^(O/i^)~^/^ Mq are the first to collapse, while smaller scales 
are of no cosmological relevance. Thus, in the isothermal picture the first masses to condense out 
of the primordial plasma have the size of a globular cluster. These first structures start clustering 
together via gravitational instability to form successively larger agglomerations. In other words, 
galaxies and galaxy clusters form hierarchically in a "bottom up" fashion. Qualitatively, one 
expects to see a roughly self-similar clustering pattern without the dramatic structures of the 
adiabatic models. 

Isothermal models do not have difficulties with the CMB anisotropics, since by definition 
the radiation field is uniformly distributed. In addition, the mass scale of the crucial first 
generation of structures is too small. The major problems with these scenarios is that the 
structure we see in the universe is not purely hierarchical. Galaxies have their own individual 
properties and galaxy clusters dot not just look like large galaxies. In the isothermal models this 
individuality of galaxies must be explained through some non-gravitational, probably dissipative, 
process. In addition, isothermal disturbances seem rather unnatural. Only very special physical 
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processes can lead to primordial fluctuations in the matter component leaving the tightly coupled 
radiation field undisturbed. Particular inflationary models, where the scalar fleld responsible for 
generating the fiuctuations is not the infiaton itself, may be able to produce this special type of 
perturbations. 

5 Non-baryonic structure formation 

The difliculties faced by the purely baryonic models opened the way for alternative structure 
formation scenarios, where the dominant matter is non-baryonic. Further motivation for con- 
sidering models dominated by exotic matter species came from a combination of observational 
data and theoretical prejudice. Observations of the light element abundances, in particular, re- 
quire Ob/i^ < 0.015 to comply with primordial nucleosynthesis calculations. On the other hand, 
dynamical considerations seem to imply that 17 ~ 0.2 and inflation argues for Q ~ 1. If all these 
are true, our universe must be dominated by exotic non-baryonic species to the extent that the 
baryons are only a small fraction of the total matter. 

5.1 Non-baryonic cosmic relics 

An ongoing problem of the dark matter scenarios is that we still do not exactly know neither the 
nature nor the masses of the particles that make up the collisionless component of our universe. 
High energy physics theories, however, provide us with a whole "zoo" of candidates which are 
known as "cosmic relics" or "relic WIMPs". Typically, we distinguish between "thermal" and 
"non-thermal" relics. The former are kept in thermal equilibrium with the rest of the universe 
until they decouple. A characteristic example of this relic type is the massless neutrino. Non- 
thermal relics on the other hand, such as axions, magnetic monopoles and cosmic strings, have 
been out of equilibrium throughout their lifetime. Thermal relics are further subdivided into 
"hot" relics, which are still relativistic when they decouple, and "cold" relics which go non- 
relativistic before decoupling. A typical hot thermal relic is a light neutrino with ~ 10 eV. 
The best motivated cold relic is the lightest supersymmetric partner of the standard model 
particles, with a predicted ^wimp ~ 1- Note that thermal relics with masses around 1 keV are 
usually referred to as "warm dark matter" . Right-handed neutrinos, axinos and gravitinos have 
all been suggested as potential warm relic candidates. 

5.2 Evolution of the Jeans mass 

If the dark component is made up of weakly interacting species, the particles do not feel each 
other's presence via collisions. Each particle moves along a spacetime geodesic, while pertur- 
bations modify these geodesic orbits. One can study the response of the of the dark matter 
component by invoking an effective pressure and treating the collisionless species as an ideal 
fluid. The associated Jeans length is obtained similarly to the baryonic one. When dealing with 
a collisionless species, however, one needs to replace Eqs. (||), @ with the Liouville equation 
(see Coles & Lucchin (1995) §). Then, 

Aj = fDMW7;-, (77) 
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where now vdm is the velocity dispersion of the dark matter component. The corresponding 
Jeans mass is 

Mj = |7rpDM(Y) , (78) 
where pdm is the density of the non-baryonic matter. 
5.2.1 Hot thermal relics 

They decouple while they are still relativistic, that is tdec < ^nr, where we assume that ^ ^eq- 
Throughout the relativistic regime vbm ~ 1, /o — P7 oc and /Odm implying that 

Aj'^^ oc and Mj^^ oc a^. Once the species have become non-relativistic and until matter- 
radiation equality, vbm oc a~^, /j ~ oc and p^M oc a~^. Recall that the particles have 
already decoupled, which means that Tdm 7^ T^- Consequently, Aj'^^ oc a and M^^^ = constant. 
After equipartition v^m oc a"^ and p ~ Pdm oc a~^, which translates into Aj*^^ oc a^^^ and 



J 

oc a"'^/'^. Overall, the Jeans mass of hot thermal relics evolves as 

2 



a 



Mj^^ oc i constant z^j- > z > Zcq , (79) 

Clearly, Afj*^^ reaches its maximum at Znr. In fact, the highest possible value corresponds to 
particles with z^r = Zeq such as neutrinos with m,^ ~ 10 eV. In this case (Mj'^^)max — 3.5 x 
10^5(17^/i2)-2 Mq (see Coles & Lucchin (1995) §). For a typical hot thermal relic (Aff ^)inax ~ 

5.2.2 Cold thermal relics 

Cold thermal relics decouple when they are already non-relativistic (i.e. t^r < tdec < ^eq)- Thus, 

for t < tjir we have vdm ~ 1, Pdm oc and p ~ cx a~^, implying that Aj'^'* oc and 

Mj^^ (X a?. In the interval between tnr and tdoc the key variables evolve as vdm oc a^^^"^ (recall 

that Tdm — until tdcc)) Pdm oc and p ~ oc a~^. As a result, Aj"^^ oc a'^/^ and 

M^^^ oc a^/"^. After the particles have decoupled Tdm 7^ T^, which means that vdm cc a^^ . At 

the same time pdm oc and p p^ cc a~^, ensuring that Aj'^^ oc a and m]'^-* = constant. After 

equilibrium vdm oc and p oc pdm oc a~^, implying that Aj'^^ oc a^^'^ and M^f^ oc a~^/^ In 
short, the Jeans mass of cold thermal relics evolves as 



^ 0? z> Zn 



Uf^ oc { 



^ Znr ^ Z ^ ^IfJec ; /'8n"\ 

constant Z^cc > Z > Zeq , 



Accordingly, the maximum value for Mj corresponds to species with tjec = ^eq- In other 
words, the sooner the particles cease being relativistic and decouple the smaller the associated 
maximum Jeans mass. Typically (Mj'^'')max ^ lO"*^^ M©. 
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5.3 Evolution of the Hubble mass 



By definition 



Mj^ = IvrpDM ( , (81) 



where ;9dm is tlie energy density of the cohisionless species. For t < t^r we have poM oc a~^, 
Xh oc t (X a^. During the interval tnr < * < ^eq we have pdm oc and Xh oc t oc a^. Finahy, 
when teq < t, /9dm oc and oc t oc a^/^. Overall the Hubble mass of the dark matter 
component evolves as 

Mh oc ^ Zjij. > Z > Zeq , (82) 



7,3/2 



^eq ^ ^ ■ 



Following definitions (77), ( |78D and (pD, one can easily verify that the Jeans mass and the 
Hubble mass are effectively identical as long as the relic species are relativistic. 

5.4 Dissipative effects 

The ideal fluid approximation for collisionless species holds on sufficiently large scales only. On 
small scales, the free geodesic motion of the particles will wipe out any structure. This process 
is known as "Landau damping" or "free streaming" . A proper calculation of the damping scale 
associated to free streaming requires integrating the collisionless Boltzmann equation. Here, 
we will only obtain an estimate of the effect in the case of hot thermal relics. To begin with, 
consider the coordinate (comoving) distance travelled by a free streaming particle 

.Fs=r^di, (83) 

Jo a 

where £fs = ixys is the corresponding physical (i.e. proper) distance. Clearly, perturbations 
in the dark matter component on scales smaller than ^^ps will be wiped out by free streaming. 
Integrating the above we find, for t < tnr when t;DM ~ 1 and a oc 

^1/2 

,1/2 ^,1/2 

xps = — t-'/''dt = ^t'/\ (84) 

^nr J "^nr 

on using the relation a = anr(V^nr)^^^- The above suggests that ^ps = 2t by the time the 
species have ceased being relativistic. During the interval tnr < * < ^eq, when vbm oc 
(i.e. vj)M = CL-ai/o) and a oc t^/^, the integration gives 



2tnr , /"* -2 J J. ^^Ti 

xps = \-anr a at 

(^nv Jtni 



1 + In ( — 



(85) 



where 2tnr/onr is the coordinate free streaming distance at inr- The associated physical free 
streaming scale is ips = aa^FS = 2tnr/«nr[l + ln(a/anr)]«- Finally, after equipartition udm oc 
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(i.e. VBM 



r/a still) and a oc t^/^. Thus, a further integration of Eq. 

2inr 



gives 



XFS 



Cljii- 

If 



1 +ln 



1 + ln 



^cq 



^nr 



+ On 



-2 



dt 



+ 



3fflnr^eq 



On 



1 + In ^ 



cq 



+ 



3t. 



eq 



■^eq 



oq 
1/2 



1/3 



on using the relation teqOnr = ^nrOeq- Thus, the total physical fee-streaming scale is 



^FS 



1 + In ( 



+ 



3t 



eq 



1 _ ( _23. 

a 



1/2' 



a . 



(86) 



(87) 



At late times, when a 3> fleqi the above approaches its maximum value 



^FS ^ (^Fs)max = 



5 + 21n( ^ 



where the scale factor has been normalized so that a = 1 at present. To obtain numerical 
estimates we need to identify the epoch the species become non-relativistic. Assuming that the 
transition takes place when T ~ mDM/3 we find (see Padmanabhan (1993) P]) 



(^FS) 



1^dm/i')'^^ Mpc, 



(89) 



where ttt-dm is the mass of the collisionless particles in units of 1 keV. Accordingly, the minimum 
scale that survives collisionless dissipation depends crucially on the mass of the dark matter 
species. For neutrinos, with m^^ ~ 30 eV we find (^Fs)max — 28 Mpc and a corresponding 
mass scale (MFs)max ~ 10^^ Mq. Note that a more accurate treatment, using the Boltzmann 
equation, gives ^fs — 40 Mpc. For a much heavier candidate, say rriDM — 1 keV, we find 
(^Fs)max ~ 0.5 Mpc and (MFs)max ~ 10^ Mq. In general, the lighter the dark matter species 
less power survives on small scales. 

Cold thermal relics and non-thermal relics have very small dispersion velocities. As a result, 
the maximum values of the Jeans mass and of the free streaming mass are very law. In this 
case, perturbations on all scales of cosmological interest grow unimpeded by damping processes, 
although they suffer stagnation due to the Meszaros effect until matter-radiation equality. After 
recombination the potential wells of the collisionless species can boost the growth of perturba- 



tions in the baryonic component on scales of the order of (M 



10^ Mp,. 



5.5 Scenarios, successes and shortcomings 

Historically, there have been two major non-baryonic structure formation scenarios. The "Hot 
Dark Matter" (HDM) models, where the dominant collisionless matter is in the form of hot 
thermal relics, and the "Cold Dark Matter" (CDM) models in which the baryonic component is 
either a cold thermal relic or a non-thermal species. It should be emphasized that, at present, 
pure HDM models are not considered viable and that the simplest CDM models are struggling 
to survive. Here, we will outline the key features of these two old adversaries and present some 
of the current alternatives. 
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Figure 2: Typical evolution of a perturbed scale Af = M{X) in HDM scenarios. During stage A 
(t < tent < teq) the mode grows as Sum t. Throughout stage B (tent < t < teq) the perturbation is 
frozen in due to the Meszaros effect. Finally, in stage C (tree < t) the mode becomes unstable again and 
grows as 6h esc t^/-^ Note that fluctuations with size smaller than Mps are wiped out by neutrino free 
streaming. 

5.5.1 Hot dark matter models 

Typical HDM scenarios involve thermal relics with z^j- ~ Zeq and the best motivated candidate 
is a light neutrino species with rrii, ~ 10 eV. The key feature of the perturbation spectrum 
is the cutoff at ^fs — 40 Mpc due to the neutrino free streaming. Because of this, the first 
structures to form have sizes of approximately lO'^^ Mq which corresponds to a supercluster of 
galaxies (see Figure 2). Moreover, because the scale is very large collapse must have occurred 
at relatively recent times (i.e. at z < 3). Thus, in a universe dominated by hot thermal relics, 
structure formation proceeds in a top-down fashion similar to the adiabatic baryonic models. 
Perturbations on scales as large as 10^^ Mq go no linear in a highly non-spherical way. As a 
result, they collapse to one dimensional objects resembling the pancake-like formations of the 
baryonic models (Zeldovich (1970) ||l^). Once the pancake forms and goes non- linear in one 
of its dimensions, the baryons within start colliding with each other and dissipate their energy. 
Thereby, the baryonic component fragments and condenses into smaller galaxy-sized objects. 
The neutrinos, however, do not collide neither with each other nor with the baryons. They are 
therefore unable to dissipate their energy and collapse into more tightly bound objects. They 
remain less condensed forming what one might call a "neutrino halo" around the baryons. 

Several groups have numerically simulated structure formation in a neutrino dominated 
universe (e.g. Centrella & Mellott (1982); White et al (1983) |l|]). In ah simulations one 
notices a cell-like structure, which reflects the damping scale due to neutrino free streaming. 
The structure of these simulations on scales larger than 10 Mpc is qualitatively similar to the 
voids and filaments seen in some of the redshift surveys. However, the models have problems 
reproducing the small scale clustering properties of galaxies. In particular, the HDM simulations 
can agree with the observed galaxy-galaxy correlation function only if the epoch of pancaking 
takes place at z ~ 1 or less. This seems too late to account for the existence of galaxies 
with redshift z > 1 and of quasars with z ~ 5. A possible way out is if the density field of 
the universe is traced by galaxy clusters rather than by galaxies themselves. In this case the 
mismatch between the predicted and the observed amount of small-scale structure is alleviated. 
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An additional problem of the HDM scenarios is that the baryons are shock heated as they fall 
into the pancakes. This might have raised the temperature to levels that could have prevented 
the baryonic matter from condensing or led to excessive x-ray production. 



5.5.2 Cold dark matter models 

Cold dark matter (CDM) candidates are cold thermal relics and non-thermal relics with ^ 
Zeq- For such particles the maximum damping scale is too small (<C 1 Mpc) to be of any cosmo- 
logical relevance. In the CDM models the important feature is the weak growth experienced by 
perturbations between horizon crossing and equipartition. It means that the density contrast 
increases as we move to smaller scales, or that the perturbation spectrum has more small-scale 
power. Note that the shape of the HDM spectrum in not cosmologically important, as the min- 
imum scales to survive free streaming are way too large. Thus in standard CDM scenarios the 
first objects to break away from the background expansion have sub-galactic sizes (< 10® Mq). 



These structures virialize through violent relaxation (Lynden-Bell (1967); Shu (1978) |17|) into 
gravitationally bound configurations that resemble galactic halos. At the same time the baryons 
can dissipate their energy and condense further into the cores of these objects. As larger and 
larger scales go non-linear, bigger structures form through tidal interactions and mergers. Hence, 
according to the CDM scenario structure forms in a bottom-up fashion analogous to that of the 
isothermal baryonic models. Note that the ability of baryons to dissipate allows objects of as- 
trophysical size to condense out as individual and discrete entities. The predicted properties of 
the different kinds of galaxies that form this way appear in good agreement with observations 
(Blumenthal et al (1984) [|l9|). 

The simplest CDM scenario, now known as the "Standard CDM" (SCDM) model, has critical 
density and contains only cold collisionless species (Peebles (1982); Blumenthal et al (1984); 
Davis et al (1985, 1992) |l^). Major problems with these scenarios are their failure of the 
simulations to reproduce the observed galaxy-galaxy correlation function and the overproduction 
of clusters when the modes are normalized with the COBE results. One can attempt to save 
the SCDM model by removing the excess small-scale power from the perturbation spectrum, 
an approach known as tilting of the spectrum, (Vittorio et al (1988); Bond (1992); Liddle et al 
(1992) ll^). This strategy, however, usually has disastrous effects on the acoustic peaks of the 



CMB. At the moment, the only way out for the SCDM model is a baryon density substantially 



higher than the one predicted by standard nucleosynthesis (White et al (1995b, 1996) [21|). 



5.5.3 Alternative options 

In an effort to salvage critical density, researchers have abandoned the CDM hypothesis in favor 
of a mixture of cold and hot dark matter particles (Bonommetto & Valdarnini (1984); Fang et 
al (1984); Shafi & Stecker (1984); Holtzman (1989); Schaefer et al (1989) ||). Hybrid "Cold- 
Hot Dark Matter" (CHDM) scenarios take advantage of the free streaming properties of their 
hot collisionless component to reduce the small-scale power in the perturbation spectrum. The 
original version had rJcDM = 0.7 and f^HDM = 0.3 (Davis et al (1992); Schaefer & Shafi (1992, 



1994); Taylor & Rowan-Robinson (1992); Klypin et al (1994) |2|]), but the free streaming effect 
was too strong. Currently, the preferred values for the hot dark matter contribution lie near 
^^HDM = 0.2 (Klypin et at (1995); Pogosyan & Starobinski (1995); Liddle et al (1996) H). Note 
that CHDM models require a rather uncomfortable fine-tuning to produce two particle species 
with similar cosmological densities but very different masses. 
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If one wants to retain the CDM hypothesis, namely to assume cold collisionless species only, 
the simple strategy is to reduce the matter density. This shifts matter-radiation equality to a 
later epoch, lengthening the period in which the small-scale modes have their growth suppressed. 
In view of the very different normalization of low-density universes with the COBE data, models 
with J7 ~ 0.3 ought to be easily distinguishable from those with critical density. This is not the 
case, however, since most other observations also require analogous renormalizations. 

The natural way of keeping the low-density models compatible with standard inflation is to 
introduce a cosmological constant, thus obtaining the ACDM scenarios (Peebles (1984); Turner 
et al (1984); Efstathiou et at (1990); Ostriker and Steinhardt (1995) H). Note that, at present, 
a non-zero cosmological constant is favoured by the type-la suprenovae measurements (Perlmut- 
ter et al, (1998); Schmidt et al (1998) ||2^). The presence of a non-zero A means that one can 
capitalize on the freedom to vary the rest of the parameters. At the moment it appears that 
~ 0.3 is the allowed lower limit on the density, with the optimum value near = 0.3. Note 
that one can achieve extra freedom by assuming that the cosmological constant is actually a 
decaying function of time, which is usually referred to as "quintessence" (Coble et al (1997); 
Turner & White (1997); Caldweh et al (1998) [|^]) 

An alternative way of shifting the time of equipartition is by adding extra massless species 
in the model (Bardeen et al (1987); Bond and Efstathiou (1991) |]28|| ). Recall that the present 
density in relativistic particles is not known since the neutrino background remains undetectable. 
Note, however, that the presence of relativistic species is strongly constrained by nucleosynthesis, 
which means that the extra relativistic energy must be generated later on (Dodelson et al (1994); 
McNally k Peacock (1995); White et al (1995a) fH). 

In all of the aforementioned models the perturbations are of the adiabatic type. Scenarios 
based on pure isocurvature fluctuations do not seem viable. However, an admixture of adia- 
batic and isocurvature modes has not been excluded. In fact, the COBE data do not seem to 
discriminate against an isocurvature component (Stompor et al (1996) pO|). 

Topological defects, such as cosmic strings, monopoles domain walls and textures, provide 
a very different alternative to structure formation (Vilenkin &; Shellard (1994); Hindmarsh & 
Kibble (1995) Q). They have much more predictive power than inflation, but they are also 
more technically demanding. Moreover, current observations appear unfavorable to them (Allen 
et al (1997); Pen et al (1997) |3l|). 

6 Discussion 

The question of how the observed large-scale structure of the universe developed and how galax- 
ies were formed has been one of the outstanding problems in modern cosmology. Looking back 
into the past hundred years one sees three decisive moments in the pursuit of the answer. The 
first milestone was the formulation of general relativity which provided researchers with the 
theoretical tool to probe the large-scale properties of the universe. Bubble's observations mani- 
festing the expansion of the universe was the second decisive moment, as it forced cosmologists 
to break away from the then prevailing concepts of a static and never-changing cosmos. Finally, 
the discovery of the Cosmic Microwave Background radiation by Penzias and Wilson estab- 
lished the Hot Big Bang theory, an idea advocated several years earlier primarily by Gamov. At 
last, cosmologists had a definite model within which they could tackle the structure formation 
question. Pending future observations, one could argue that the inflationary paradigm is one 
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additional milestone in our effort to understand the workings of our universe. The recent su- 
pernovac measurements, which suggest an accelerated universal expansion, could prove another 
very decisive moment. Time will show whether they actually are. 

Despite the problems and the uncertainties, cosmologists now believe that all the structure 
that we observe around us today originated from minute perturbations in a cosmic fluid that 
was smooth to the accuracy of one part in ten thousand at the time of recombination. Such tiny 
irregularities could have been triggered by quantum fluctuations that were stretched out during 
the inflationary expansion or by topological defects such as cosmic strings for example. Given the 
current observational status, inflation appears to be the strongest candidate. In these notes we 
have set aside the question of the origin of the primeval fluctuations. Our discussion focussed on 
the linear evolution of these minute irregularities once the universe entered the post-inflationary 
Hot Big Bang era, and on the physical processes that could have affected them. When studying 
density perturbations, one sooner or later encounters the reality that baryon inhomogeneities 
cannot actually grow before recombination. This fact, together with the extreme smoothness of 
the CMB temperature, implies that baryonic fluctuations simply do not have enough time to 
produce the observed structure. When one adds to that the nucleosythesis constraints and the 
strong theoretical and observational bias for spatial flatness, the once popular baryon-dominated 
picture of the universe seems unsustainable. So, if the baryons are not the dominant form of 
matter in our universe, then what is? The answer to that might lie in high energy physics 
theories. Theoretical physics provides a whole zoo of supersymmetric, dark matter species that 
could bring $7 close to unity and also "assist" structure formation. The attractive feature of 
collisionless matter is that perturbations in its density start growing earlier than those in the 
baryonic component. Thus, as soon as the baryons decouple from radiation, they undergo a 
period of fast growth as they fall into the potential wells of the collisionless species. This can 
improve the final picture but unfortunately does not solve all the problems. The key obstacle 
being that a single collisionless species does not seem capable of fitting all the data. At best, one 
needs the presence of one cold and one hot dark matter species in a combination that will make 
the most of their advantages while minimizing their shortcomings. The recent supernovae results, 
suggesting that our universe might be dominated by some sort of dark energy have added extra 
flavor to the whole situation. The presence of a cosmological constant, or quintessence, means 
that researchers have extra freedom when dealing with crucial cosmological parameters, such as 
the age of the universe for example. On the other hand, however, the idea that the majority 
of the matter in our universe is in the form of some unknown exotic species brings back some 
rather embarrassing memories from our relatively recent past. A lot could be decided within the 
next few years as we expect an influx of high quality data. The "Boomerang" and "Maxima" 
observations have added valuable information which seems to favor the inflation based models. 
The near future "Probe" and "Planck" satellite missions also promise high precision data. For 
some researchers structure formation is a story that is fast reaching its conclusion. The future 
will show if they are right or just too hasty. 
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